Steam reforming on transition-metal carbides from density-functional theory 
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A screening study of the steam reforming reaction (CH4 + H2O — > CO + 3H2) on early transition- 
metal carbides (TMC's) is performed by means of density-functional theory calculations. The set 
of considered surfaces includes the a-Mo2C(100) surfaces, the low-index (111) and (100) surfaces of 
TiC, VC, and 5-MoC, and the oxygenated a-Mo 2 C(100) and TMC(lll) surfaces. It is found that 
carbides provide a wide spectrum of reactivities towards the steam reforming reaction, from too 
reactive via suitable to too inert. The reactivity is discussed in terms of the electronic structure of 
the clean surfaces. Two surfaces, the <5-MoC(100) and the oxygen passivated a-Mo2C(100) surfaces, 
are identified as promising steam reforming catalysts. These findings suggest that carbides provide 
a playground for reactivity tuning, comparable to the one for pure metals. 



I. INTRODUCTION 

Steam reforming is an important industrial process, 
where natural gas (CH4) is converted into synthesis gas 
(CO and H2) according to the overall reaction 



CH 4 + H 2 ^ CO + 3H 2 



(1) 



The synthesis gas is subsequently transformed into more 
valuable chemicals, such as ammonia, methanol, and 
diesel. Since steam reforming acts as a source of hy- 
drogen, it is also potentially important for any emerging 
hydrogen economy. For a detailed review of the steam re- 
forming process and its applications the reader can con- 
sult Ref. Q. 

Commercially the steam reforming reaction is con- 
ducted over a Ni-based catalyst due to the relatively 
low cost and good activity of nickel. This reaction has 
been studied in detail on the close-packed Ni(lll) sur- 
face and on the stepped Ni(211) surface by means of 
density-functional theory^ The major technological chal- 
lenge for Ni catalysts is the formation of carbonaceous 
deposits, termed coke, that lead to catalyst deactivation. 
In Refs. it was established that the step edges on 
Ni surfaces act as growth centers for graphite. Other 
transition metals (TM), such as Ru, Rh, Pd, Ir, and 
Pt also show high activity and selectivity towards steam 
reforming 4 and have a high resistance against carbon for- 
mation. However, these materials are scarce in nature 
and expensive. Therefore new materials that are resis- 
tant to carbon formation are needed. 

Transition-metal carbides (TMC's) have gained quite 
some attention since Levy and Boudart reported that 
they have "platinum- like behavior" for certain reactions. 5 
The starting material for the production of carbides is 
cheap and abundant and therefore it has been suggested 
that they can replace the noble metals in catalysis. The 
main problem with the carbides has been to produce ma- 
terials with sufficiently high surface area for them to be 
interesting for catalytic applications. This problem has 
been overcome and several studies report on carbides 
with surface areas as high as 200 m /g i ' ' ' 

In Refs. [13 and [n]> it was shown that carbides of 
molybdenum and tungsten are stable and extremely ac- 



tive catalysts for not only the steam reforming but also 
the dry reforming and the partial oxidation of methane at 
elevated pressure. Moreover, no macroscopic carbon was 
deposition on the catalysts during the catalytic reactions. 
The relative activity of a number of steam reforming cat- 
alysts was established as: Ru > Rh > Ir « Mo 2 C > WC> 
Pd > Pt. 

Several density-functional theory studies concerning 
adsorption on TMC's have been conducted, for example, 
atomic adsorption,i2ii2iIli5iI£LiLI&i2i2a 2 adsorption^ 
CO adsorption, 2 ^ 2 ^ jsjjj^ = l - 3) adsorption, 1 ^ 
methanol adsorption, 25 CH K (x = — 3) and C 2 H4 
adsorption*^ In addition, the water-gas-shift reaction 
has been studied in Refs. [^IHIM 

To our knowledge, a theoretical understanding of the 
steam reforming reaction on TMC's is missing in the lit- 
erature. The aim of this paper is to investigate how the 
TMC's are suited for the steam reforming reaction by 
performing a set of density-functional theory (DFT) cal- 
culations on different surfaces of a-Mo 2 C, TiC, VC, and 
S-MoC in order to catch the trends between TMC's com- 
posed of different TM atoms and of different structures. 
The influence of the surfaces being covered with oxygen 
is also investigated. 

The outline of this paper is as follows. To begin with 
the systems to be investigated are defined in Sec.|TTJ This 
is followed by Sec. lIIIi which describes the calculation de- 
tails. Our results are presented in Sec. IIVI and discussed 
in Section [V] The paper ends with a conclusion in Sec- 
tion ED 



II. INVESTIGATED TMC SYSTEMS 

The Mo-C system can occur in several crystalline 
forms, see for example the phase diagrams in Refs. HO 
and [3l|. There are two types of Mo 2 C phases, one or- 
thorhombic (a) and one hexagonal (/?), which is a high 
temperature phase. In this study, we focus on the a 
phase. For this phase the Mo atoms are only slightly dis- 
torted from an hep arrangement and the carbon atoms 
occupy one half of the octahedral interstitial sites. In 
the present molybdenum carbide literature there seems 
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FIG. 1: Top view of the relaxed (a) Mo-terminated 5-MoC(lll), (b) C-terminated <5-MoC(lll), (c) <5-MoC(100), (d) Mo- 
terminated q-Mo2C(100), (e) C-terminated a-Mo2C(100) surface. The blue and black balls represent the Mo and C atoms, 
respectively. The black dashed lines indicate the unit cell. 



to be a confusion when it comes to the attribution of the 
symbols to the different phases, especially regarding the 
a and j3 phases. This is probably because the a-Mo2C 
phase was originally assigned a hexagonal structure but 
refined experiments identified the structure to be of or- 
thorhombic type 3 -^, with lattice parameters a — 4.729 A, 
b = 6.028 A and c = 5.197 AM 

The close-packed surface of CK-M02C is the (100) sur- 
face. It consists of alternating Mo and C layers and 
can be either Mo- or C-terminated. Experimental stud- 
ies of a-Mo2C(100) show that the surface termination 
and structure depend strongly on the cleaning procedure. 
Both the Mo- and the C-terminated surfaces have been 
observed i 34 ' 35 Therefore we consider both the TM- and 
the C-terminated carbide surfaces. 

The TiC, VC and S-MoC adopt the NaCl structure. 
The (100) surface of these carbides consists of layers of 
equal amounts of TM and C atoms. The (111) surface, 
however, is polar and consists of alternating TM and C 
layers, similar to the a-Mo2C(100) surfaces. One im- 
portant difference is that the C-terminated a-Mo 2 C(100) 
surface has half the amount of C atoms as compared to 
a C-terminated TMC(lll) surface. 

In this paper, we report the results from DFT cal- 
culations on (i) the TM- and C-terminated surfaces of 
a-Mo 2 C(100), TiC(lll), VC(lll), and 5-MoC(lll); (ii) 
the O-covered surfaces in (i); and (iii) the (100) surfaces 
of TiC, VC and S-MoC. 

III. THEORETICAL DETAILS 

All the DFT calculations in this paper are performed 
with the computer code DACAPO 3 -^, which uses plane- 
waves and ultra-soft pseudopotentials. For the optimiza- 
tion of the bulk lattice parameters of the considered 
TMC's a plane-wave cutoff of 400 eV, a 8 x 8 x 8 k- 
point sampling, and the RPBE exchange-correlation (xc) 
functional 3 -^ are used. The bulk structure of a-Mo2C was 
obtained by minimizing the total energy of the unit cell 
with respect to the length of the lattice vectors using a 
Newton-Raphson scheme and allowing for complete in- 
ternal relaxation at each step. The calculated lattice pa- 
rameters a = 4.825 A, b = 6.162 A, and c = 5.304 A 



differ by ~ 2% from the experimental values^ 3 - 

For TiC, VC and S-MoC, the equilibrium bulk struc- 
ture was obtained utilizing the Murnaghan equation of 
stated The calculated lattice constants are 4.336 A, 
4.171 A, and 4.444 A for TiC, VC, and <5-MoC, respec- 
tively. A quantitative agreement is found between the 
theoretical and experimental lattice constants ! 39 ' 40 

The surfaces are modeled using the supercell approach, 
with the DFT lattice parameters given above. A 400 eV 
energy cutoff for the plane-wave expansion of the one- 
electron orbitals, a 4 x 4 x 1 fc-point grid, and the RPBE 
xc functional are used. Each TMC(lll) and TMC(100) 
surface is represented by a slab with the geometry defined 
in Fig. Q] and a thickness of four atomic layers. Repeated 
slabs are separated by at least a 10.8 A thick vacuum 
region. The Mo- and C-terminated a-Mo2C(100) slabs 
have a surface geometry as defined in Fig. [T] and consist 
of four layers Mo and four layers C separated by 10.6 A 
vacuum. During structural relaxations, the two bottom 
atomic layers are constrained to the corresponding bulk 
geometry. 

Figure Q] shows the relaxed structures of the Mo- and 
C-terminated J-MoC(lll) surfaces, the J-MoC(lOO) sur- 
face, and the Mo- and C-terminated a-Mo2C(100) sur- 
faces. These relaxed surfaces and the corresponding VC 
and TiC surfaces are used as starting points for the mod- 
eling of the steam reforming reaction. 

The investigated steam reforming reaction is assumed 
to have the following nine elementary steps: 



CH 4 + 2* 




CH*. 


+ H* 


(2) 


ch*,+* 
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+ H* 
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CO* 
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2H* 




H 2 + 


2*, 


(10) 



where * (2*) denotes one (two) empty surface site and X* 
denotes the adsorbate X bonded to this site. Adsorption 
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of all the intermediates (CH 3 , CH 2 , CH, C, OH, O, CO, 
and H) in Eqs. (p?|)- (fITj)) are considered. In all adsorption 
calculations, the adsorbate coverage is one-quarter of a 
monolayer (ML). Adsorption is allowed on one of the two 
slab surfaces at the time. Adsorbates together with the 
surfaces are allowed to relax in all directions until the 
forces are less than 0.05 eV/A. The energies in the cal- 
culated potential energy diagrams are expressed relative 
to the clean surface (O-covered surfaces if the surface is 
oxygenated), a CH4 molecule, and a H 2 molecule in the 
gas phase. 

The steam reforming reaction on the oxygenated TMC 
surfaces is modeled by adsorbing the reaction interme- 
diates on a 1ML O-covered TMC surface. According 
to DFT calculations in Ref. [28|, the O covering of the 
Mo- and C-terminated Mo2C(100) surfaces is energeti- 
cally favorable up to 1ML. Our calculations give the same 
structure for the 1ML O-covered Mo- and C-terminated 
Mo 2 C(100) surfaces as in Ref. HI 



IV. RESULTS 

This Section focuses on the results obtained from the 
screening calculations study. We provide the calculated 
potential energy diagrams for the steam reforming re- 
action on all the considered TMC surfaces. These are 
compared with existing results on Ni(lll) and Ni(211) ? . 
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A. TMC(lll) and a-Mo 2 C(100) surfaces 

The calculated potential energy diagram for the steam 
reforming reaction on the TM- and C-terminated TMC's 
is shown in Fig. [2] First of all we notice that the span of 
the potential energies of the investigated surfaces is huge. 
Compared to the Ni surfaces^, all the TMC surfaces are 
much more reactive because they bind the intermediates 
too strongly. Therefore none of these surfaces is suitable 
for the steam reforming reaction. 

In the following we present the trends in reactivity 
that are found when changing the TM atom in the TMC 
compound. For each of the considered TMC's, the C- 
tcrminated surface is found to be more reactive than the 
TM-terminated one. This is to be expected since the 
TMC(lll) surfaces are found to be TM terminated under 
UHV conditions" All the TM-terminated TMC (111) 
surfaces have similar potential energy shape as a func- 
tion of the reaction intermediates. The same is found for 
the C-terminated group of TMC surfaces. For the TM- 
terminated TMC surfaces, the reactivity can be ordered 
as TiC(lll) > VC(lll) > a-Mo 2 C(100) > (5-MoC(lll), 
while for the C-terminated TMC surfaces we find that the 
reactivity order is TiC(lll) > VC(lll) > S-MoC(lll) > 
a-Mo 2 C(100). Hence different trends are observed for the 
different terminations. 

The considered a-Mo 2 C(100) surfaces are less reac- 
tive than the C-terminated <5-MoC(lll) surface but more 



FIG. 2: Calculated potential energy diagram for the steam 
reforming reaction on the TM-terminated (circles) and C- 
terminated (squares) TMC(lll) and a-Mo2C(100) surfaces. 
For comparison the Ni(lll) and Ni(211) data, adapted from 
Ref. |2|, are given. 



reactive than the Mo-terminated <5-MoC(lll). For a 
given carbide, the smallest difference in reactivity be- 
tween the TM- and the C-terminated surfaces is found 
for a-Mo 2 C(100). This can be due to the fact that the 
C-terminated a-Mo 2 C(100) surface has 50% less C than 
the C-terminated TMC(lll) surfaces. 

The largest minimum in the potential energy surface is 
found upon adsorption of O-carrying intermediates, that 
is, a strong O-surface bond is formed. This implies that 
the TMC(lll) surfaces are easily oxidized. Experimen- 
tal studies show that the 2 adsorption on TiC(lll) 43 
and a-Mo 2 C(100) 3 — are dissociative. Therefore we in- 
vestigate how the steam reforming reaction is influenced 
when the TMC(lll) and the a-Mo 2 C(100) surfaces are 
oxygenated. 



B. O-covered TMC(lll) and O-covered 
a-Mo 2 C(100) surfaces 

Figure[3]shows the calculated potential energy diagram 
for the steam reforming reaction on the 1ML O-covered 
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FIG. 3: Calculated potential energy diagram for the 1ML O- 
covered TM-terminated (circles) and C-terminated (squares) 
TMC(lll) and a-Mo 2 C(100) surfaces. Also given are the 
Ni(lll) and Ni(211) data, adapted from Ref.0 



TM- and C-terminated TMC surfaces. We find that the 
O-covered C-terminated a-Mo2C(100) has a potential en- 
ergy profile similar to that reported for the Ni(lll) sur- 
face, which is the preferred Ni facet for the steam reform- 
ing reaction according to the DFT calculations in Ref. @- 
The other O-covered surfaces are unsuitable as catalysts 
for the steam reforming either because several of the reac- 
tion steps are significantly uphill energetically, or because 
the surface gets poisoned by some of the intermediates. 

Some of the trends that are found for this class of sur- 
faces are here described. Each of the O-covered TMC sur- 
faces is less reactive than its corresponding non-oxidized 
surface. For a given O-covered TMC surface the TM- 
terminated surface has higher reactivity than the C- 
terminated one, in contrast to the finding on the non- 
oxidized surfaces. For the O-covered TM-terminated 
TMC surfaces, the reactivity can be ordered as follows 
a-Mo 2 C(100) > <S-MoC(lll) > VC(lll) > TiC(lll), 
while for the O-covered C-terminated TMC surfaces we 
find that a-Mo 2 C(100) > TiC(lll) > VC(lll) > 8- 
MoC(lll). The order in reactivity between the differ- 
ent surfaces is thus different for the TM- and for the 
C-terminated O-covered surfaces. 
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FIG. 4: Calculated potential energy diagram for TMC(IOO). 
The Ni(lll) and Ni(211) data given here are adapted from 
Ref. 0. 



C. TMC (100) surfaces 

The calculated potential energy diagram for the steam 
reforming reaction on the TMC(IOO) surfaces is shown in 
Fig-IU The i5-MoC(100) surface is found to have a poten- 
tial energy curve similar to the one of the Ni(lll) surface. 
We find the same shape of the potential energy curves for 
all the TMC(IOO) surfaces. A local minimum is found for 
all the TMC(IOO) surfaces upon adsorption of CO. The 
reactivity of the surfaces follows the order 5-MoC(100) > 
TiC(lOO) > VC(100), that is, a non linear variation when 
moving from right to left in the periodic table. Com- 
pared to the TMC(lll) surface the TMC(IOO) surfaces 
are far less reactive. In Ref. HI it was suggested that the 
TiC(lOO) surface is a suitable catalyst for the water-gas- 
shift reaction. However, according to our calculation it is 
not the surface of choice for the steam reforming reaction 
since it is inert (see Fig. [4|. 



V. DISCUSSION 

The conducted computational screening study clearly 
shows that only two surfaces, the O-covered C- 
terminated a-Mo 2 C(100) and the <5-MoC(100), are suit- 
able for the steam reforming reaction. These two surfaces 
are found to have potential energy profiles similar to the 
one of Ni(lll), which is the preferred Ni surface for steam 
reforming according to Ref. Q. 
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The obtained potential energy diagrams show that 
there is a large variety in the reactivities of the differ- 
ent carbides (see Figs. [21 [3] and H}. However, besides 
the two candidate surfaces none of the other considered 
TMC surfaces are suitable for the steam reforming reac- 
tion since they either get self-poisoned by the intermedi- 
ates or they show a significant energetical uphill for some 
of the elementary steps in Eqs. ([2"])-Q10p. 

An understanding of the reactivity can be achieved 
by stud ying th e electronic structure of the surface. In 
Refs. EllllJUmiiS a tllorou g n analysis of the elec- 
tronic structure established that the TM-terminated 
TMC(lll) surfaces possess surface resonances (SR's). 
These were identified to be responsible for strong 
adsorbate-surface bonds. Therefore the high reactivities 
of the TM-terminated TMC (111) surfaces (see Fig. |2J) 
are attributed to the presence of these SR's. Preliminary 
electronic structure calculations show that similar reso- 
nances are present also on the layered polar a-Mo 2 C(100) 
surfaces, which would explain their high reactivities. 

The NaCl TMC(IOO) surfaces, on the other hand, do 
not posses any SR' s 17 ' 44 and show a much lower reac- 
tivity (see Fig. [3} compared to the TMC (111) surfaces. 
Since a suitable candidate, the 5-MoC(100) surface, is 
found among these surfaces we have tried to change the 
reactivity of the TMC (111) surfaces by weakening their 
SR's. 

The lower reactivities of the O-covered TMC(lll) and 
a-Mo2C(100) surfaces (see Fig. [3]), compared to the cor- 
responding clean surfaces can be understood as follows. 
Adsorption of oxygen on the TMC(lll) results in a 
quenching of the SR's. Thus, the surfaces are passivated, 
that is, their ability to form strong bonds with other 
adsorbates decreases due to the strong O-surface bond. 
This result shows that the reactivity of carbides surfaces 
can be tuned. 

An important question for catalytic applications is the 
stability of the candidate materials. Steam reforming op- 
erates over a large range of working conditions with tem- 
peratures 700 — 1100 °C and pressures 1 — 25 bar. The 
most recent phase diagram for Mo-C systems is the one 
by Velikanova et alM- It shows that the bulk CK-M02C 
phase exists for temperatures < 1440 °C, while the 8- 
MoC is a high temperature phase which exists above 
1956 °C. Based only on bulk data, a-Mo2C is preferred 
over (5-MoC. 

When it comes to the stability of a surface the situation 
can of course differ from the bulk. According to York and 
coworker s 10 ' 11 , M02C was active and selective for the sto- 
ichiometric steam reforming of methane to synthesis gas. 
The catalyst deactivated at atmospheric pressure, but 
was very stable when elevated pressures were employed, 
and no carbon deposition was observed on the catalysts. 
In those studies the surface characteristics was, however, 



not addressed. 

A surface characterization of a-Mo2C(100) has been 
performed in Ref. [U. The surface preparation consisted 
of Ar ion bombardment followed by annealing at suc- 
cessively higher temperatures. The surface was found to 
change from being Mo terminated to C terminated at an- 
nealing temperatures above 1300 K. Despite that these 
surfaces are stable at temperatures relevant for the steam 
reforming, our study shows that they are too reactive. 
Adsorption of O2 on both the Mo- and the C-terminated 
a-Mo2C(100) surfaces were found to be dissociative and 
the surfaces were covered with oxygen at 800 K. 35 This 
supports the stability of the O-covered a-Mo 2 C(100) sur- 
face. 

In addition, stable crystalline molybdenum carbide 
films prepared by chemical vapor deposition (CVD) and 
physical vapor deposition (PVD) have been found to con- 
tain either Mo 2 C and/or the (J-MoCi-^ phase4& 1 4 ML2& 
Hence it seems plausible that both the O-covered C- 
terminated a-Mo 2 C(100) and the <5-MoC(100) surfaces 
may be stable at relevant steam reforming conditions, 
however, further studies are needed. 



VI. CONCLUSION 

In this paper, a computational screening study of the 
steam reforming reaction on transition-mental carbides 
has been presented. By using a screening approach we are 
able to map out the huge reactivity space, which ranges 
from very reactive to inert, provided by the carbides. A 
small set of surfaces, consisting of the <5-MoC(100) and 
the O-covered C-terminated a-Mo2C(100), is identified 
as promising candidate for the steam reforming reaction. 
With this thermochemical data at hand we lay a basis for 
a thermodynamic and micro-kinetic treatment of the car- 
bides. This study illustrates the versatility and tunabilty 
of the TMC systems for the steam reforming reaction. 
The results can be understood in terms of the previously 
proposed model of chemisorption on TMC surfaces, in 
which surface resonances play crucial roles. However, fur- 
ther studies of the surface electronic structure are needed 
to establish what surface characteristics that are required 
of a carbide to make is a suitable catalyst for steam re- 
forming. 
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